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INTRODUCTION 


The  dynamic  modeling  part  of  an  experimental/theoretical  effort  to  develop  and 
test  pointing  and  tracking  control  techniques  for  a  mechanical  system  containing  both 
rigid  and  flexible  bodies  is  presented  in  this  report.  Recently,  much  work  was  done  to 
model  open-loop  chains  of  ridig  and  elastic  bodies  with  application  to  kinematic 
linkages,  spacecraft,  and  manipulators. 

In  this  modelling  effort,  the  work  of  Schmitz1  was  used  as  a  guide.  The  mechani¬ 
cal  system  consisted  of  a  flexible  beam  fixed  to  a  rigid  disk  and  rotating  in  the  horizontal 
plane.  The  new  contributions  that  result  from  this  work  are:  (1 )  the  inclusion  of  the  size 
(the  radius  of  the  rigid  disk)  as  well  as  the  mass  properties  of  the  rigid  moving  base  in 
deriving  the  equations  of  motion,  (2)  the  study  of  the  effect  of  the  size  of  the  moving 
rigid  base  on  the  dynamic  response  of  the  system,  and  (3)  the  development  and  com¬ 
parison  of  several  modelling  approaches  which  can  be  used  for  control  design:  a 
reduced-order  model  with  cantilever  modes,  a  reduced-order  model  with  system  modes, 
a  model  based  on  finite  element  analysis,  and  a  lumped  parameter  model. 

The  computer  software  used  for  this  work  was:  MSC  Pal  2  (finite  element 
analysis),  MATRIXx  (dynamic  analysis  and  control  design),  and  MathCAD 

(multipurpose  engineering  software).  The  book  of  Meirovitch2  was  referred  to  often 
during  the  course  of  this  work. 

Work  currently  engaged  in  is:  (1)  completion  of  model  validation,  including  study 
of  the  effects  of  Coulomb  and  viscous  damping;  (2)  development  and  testing  of  several 
control  algorithms:  PID,  LQG/LTR,  adaptive  and  Hq;  and  (3)  addition  of  an  end-point 
position  sensor  to  the  experimental  apparatus  to  allow  for  end-point  position  feedback 
control. 


Schmitz,  E.,  "Experiments  on  the  End-Point  Position  Control  of  a  Very  Flexible  One- 
Link  Manipulator,"  Ph.D.  Thesis,  Stanford  University,  Guidance  and  Control  Laboratory, 
Department  of  Aeronautics  and  Astronautics,  Stanford,  CA,  1985. 

2Meirovitch,  L.,  Analytical  Methods  in  Vibrations,  The  Macmillan  Company,  New  York, 
1967. 
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EXPERIMENTAL  APPARATUS 


An  experimental  apparatus  (fig  1)  was  built  to  demonstrate  pointing  and  tracking 
control  strategies  for  fast-moving,  lightweight,  flexible  appendages  to  rigid  bodies.  It 
consists  of  a  one  degree  of  freedom  rigid  disk,  rotating  about  a  vertical  axis,  with  a  long, 
thin,  flexible  beam  fixed  to  it  in  the  horizontal  plane.  Two  dc  torque  motors,  one  used 
for  disturbance  input  and  one  used  for  control,  cause  the  disk  to  rotate  about  the  vertical 
axis.  A  resolver  is  used  to  measure  disk  angle,  which  is  the  only  sensor  on  the  system 
at  the  present  time.  A  beam  end-point  position  sensor  is  presently  being  added  to  the 
system. 


DYNAMIC  MODELLING 

Dynamic  equations  were  developed  for  a  mechanical  system  consisting  of  a 
flexible  beam  fixed  to  a  rigid  disk.  Goals  were  to  understand  the  fundamental  charac¬ 
teristics  of  the  system  and  to  develop  a  model  of  the  system  which  can  be  used  to  test 
new  control  algorithms. 

Model  Description  and  Derivation  of  Equations  of  Motion 

Consider  a  uniform  flexible  beam  of  length  L,  fixed  to  a  rigid  disk  of  radius  R,  and 
moving  in  the  horizontal  plane,  as  shown  in  the  figure  below. 


u(t)  =  externally  applied  torque 

J  =  mass  moment  of  inertia  of  rigid  disk  about  axis  of  rotation 
R  =  radius  of  rigid  disk 
E  =  modulus  of  elasticity  of  beam  material 

I  =  area  moment  of  inertia  of  beam  cross  section  about  neutral  axis 
A  =  cross-sectional  area  of  beam 

p  =  density  of  beam  material  a  =  p  A  =  mass  per  unit  length 
L  =  length  of  beam 
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Displacement  of  any  point  P  along  the  beam’s  neutral  axis  at  a  distance  x  from  the 
disk  is  given  by  the  arbitrarily  large  disk  angle  0(t)  and  the  small  elastic  deflection  w(x,t) 
measured  from  the  line  OX.  The  assumptions  made  in  deriving  the  equations  of  motion 


t  Neglect  axial  deformations  and  axial  forces.  The  dynamic  stiffening  effect 
and  increased  natural  frequencies  that  result  from  the  axial  forces  are  assumed  negli¬ 
gible  in  this  situation. 

•  Contribution  of  the  rotational  moment  of  inertia  of  an  elementary  beam 
section  to  the  kinetic  energy  of  the  system  is  neglected.  This  is  valid  if  the  cross- 
sectional  dimensions  of  the  beam  are  small  compared  to  its  length. 

•  Effects  of  shear  deformation  are  negligible.  This  is  valid  if  the  cross- 
sectional  dimensions  of  the  beam  are  small  compared  to  its  length  and  if  the  frequency 
of  beam  oscillations  is  low. 

•  Elastic  deflection  of  the  beam  and  angular  velocity  of  the  disk  are  small. 

•  Neglect  any  energy  dissipation.  This  will  be  added  later  in  the  analysis. 

Expressions  for  the  kinetic  energy,  T,  and  the  potential  energy,  V,  of  the  mechani¬ 
cal  system  are  as  follows: 

r  -]2  r  -i2 

T=j  J  eft  e(t)  +  l  Jo  °  (R  +  x)  ■  ^e(t)  +  ^w(x,t)  dx 
r  2  i2 

V  =  1  •  Jo  E  •  I  •  ^  w(x,t)  dx  -  U(t)  •  0(t) 

Given  the  expressions  for  the  kinetic  and  potential  energy,  the  Hamilton’s  principle 
was  applied: 

8  •  J!,  (T  -  V)  dt  =  0 

The  equation  of  motion  and  boundary  conditions  that  result  are: 

Equation  of  motion: 

El-  ~ 3  w(x,t)  +  a  •  (R  +  x)  •  0(t)  +  w(x,t)  =  0 
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Boundary  conditions 

•  J  [gp  9(0  U(l)-E  I  ^  w(x,t)  E  I  ^  ■  w(x,t)J  ^0 

•  w(x,t)  =  0  at  x  =  0  (deflection  at  fixed  connection  =  0) 

— j 

•  ^  •  w(x,t)  =  0  at  x  =  0  (slope  at  fixed  connection  =  0) 

•  ^  ■  w(x,t)  =  0  at  x  =  L  (moment  at  free  end  =  0) 
d3 

•  •  w(w,t)  =  0  at  x  =  L  (shear  force  at  free  end  =  0) 

Definition  of  a  new  coordinate:  y(x,t)  =  (R  +  x)  ■  0(t)  +  w(x,t)  =  inertial  lateral 

displacement  of  the  beam  at  a  distance  x  from  the  edge  of  the  disk. 

The  equation  of  motion  and  boundary  conditions  reduce  to  the  following: 
e  i  ~,y(x,t)  +o-  £  y(x.t)  =o 

•  J  Jp-9(t)  -  U(t>  -El  ^y(x.l)  +R  El  J^y(x,t)  o=0 

•  y(x,t)  =  R  •  0(t)  at  x  =  0 

•  dx  '  y(x,t)  =  9(t) 

•  dx5  '  y(x,t)  =  0  at  x  =  L 

•  dx3  '  y(x,t)  =  0  at  x  =  L 

The  following  section  show  several  ways  to  solve  these  equations  in  a  form  which 
is  useful  for  feedback  control  systems. 


4 


Exact  Solution  and  Open>Loop  Transfer  Functions 

The  variables  y(x,s),  u(s),  and  d(s)  are  defined  as  the  Laplace  transforms  of  the 
variables  y(x,t),  u(t),  and  B(t),  respectively.  Taking  the  Laplace  transform  of  the  equa¬ 
tion  of  motion  and  boundary  conditions,  the  following  is  obtained: 

Equation  of  Motion 

+  <£TSj -  •  y(X,S)  =  o 
Boundary  Conditions: 


d4 

d'x* 


y(x,s) 


J  s2  0(s)  -  u(s)  -El 


_d2 

dx2 


y(x,s) 


+  R  •  E  I 


x  =  0 


d3 

dx3 


y(x,s) 


=  0 


x  =  o 


•  y(x,s)  =  R  0(s)  at  x  =  0 

*  dx  y(x’s)  =  e(s)  at  x  =  0 
d2 

*  dx5  y(x,s)  =  0  at  x  =  L 
d3 

•  dx3  y(x-s)  =  o  at  X  =  L 


The  following  dimensionless  ratios  are  defined: 


R 

L 


=  8 


J 


=  e 


Note  that  o  L3  =  (o  L)  •  L2  =  m  L2  where  m  is  the  total  mass  of  the  beam.  The  dimen¬ 
sionless  complex  number  X  related  to  the  Laplace  variable  s  is  defined  by: 


The  general  solution  of  the  equation  of  motion  is: 

y(x,s)  =  A  ■  sin((3  ■  x)  +  B  ■  sinh((3  ■  x)  +  C  cos((3  •  x) 
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found: 


With  this  equation  and  the  boundary  conditions  stated  above,  the  following 


r a' 

i 

o 

r  0- 

M{X)  ■ 

c 

o 

_0_ 

iL 

E  1  |3 

L-u(s)J 

D  =  R  •  0-  C 

e 

B  =  _ -  A 

P 


[-sin(X)  -  sinh(X)]  (-cos(X)  -  cosh(X)]  R  •  cosh(X)  + 

r 

[-cos(X)  -  cosh(X)]  [sin(X)  -  sinh(X)]  — +  R  •  sinh(X) 

3 

2-  5  X  -2  £  •  iL  .  LA  .  R 

P  P 


=  M(k) 


Solving  for  A,  C,  and  0: 


A_  _L  r,u(S)  ~ 

K  !  2 


_1_ 

P 


[1  +  sin(X)  •  sinh(>i)  +  cosfi)  •  cosh(^)] 


|_E  •  I  ■  (3  j  L  +  (-R)  [cos(X)  -  sinh(\)  +  sin(X)  •  cosh(X)] 


c=  _L.  ~  u(s) 

K  2 

LE  I  •  |3  j 


r  R  [1  +  cos(X)  ■  cosh(A.)  -  sin(X)  •  sinh(A.)]  -  •  • 

|  +  •  [cos(\)  sinh(^)  -  sin(\)  •  cosh(^)] 

L  P  j 


JL  .  2  u(s) 

r 
C 

E  •  I  ■  P 


0  =  —  '  -  uw  ■  [1  +  cos(X)  ■  cosh(^)] 

K  2 


9  r  3  r  ?  ? I  1 

K  =  -=-  •  [1  +  cos(X)  •  cosh(?i)3  ■  e  ■  X  +  sin(X)  •  cosh(X)  1 1  +  X  •  5  J  •  •  •  ' 

P  L  j 


!  p  p 

+  2  ■  X  ■  5  •  sin(?i)  ■  sinh(A.)  +  cos(^)  •  sinh(\)  •  '  k  -5-1 

L 


J 


The  following  open-loop  transfer  functions  are  of  interest: 


0(S) 

u(s) 


Disk  angle  sensor  colocated  with  the  actuator. 


yM ) 

e(s) 


Tip  position  sensor  (now  being  added  to  experimental 
apparatus)  non-colocated,  i.e. ,  separated  from  the  actuator 
by  the  elastic  beam. 


B(s)  N0(X)  .  y(L,s)  _  Ntip(X) 

u(s)  =  D(Xy  and  u(sr  ~  'W 


_  *  ^  2  2 
D(X)  =  [1  +  cos(X)  •  cosh(X)]  •  e  •  X  +  sin(A.)  cosh(X)  •  1  +  X  8  •  ■  ■  ’ 

L+  2  ■  X  •  5  •  sin(X)  ■  sinh(A.)  +  cos{X)  •  sinh(X)  ■  £  •  5^  -  1  | 

r  2 

N0(A.)  =  ,  -2  — 2  •  [1  +  cos(X)  •  cosh(^)] 

Le  i  xj 

r  2 

Ntip(/V)  =  - 2  (L  •  [sin(X)  +  sinh(X)]  +  R  ■  X  •  [cos(A.)  +  cosh(A.)]) 

3 

i_E  •  i  X_ 

The  roots  of  the  function  D(^)  are  the  poles  of  the  transfer  functions,  also  called 
the  natural  vibration  frequencies  of  the  rigid  disk/flexible  beam  system.  As  stated  in  the 
initial  assumptions,  this  equation  is  valid  only  to  predict  the  first  few  vibration  frequen¬ 
cies  of  the  system. 

D(X),  N0(X),  and  Ntip(X)  are  expressed  in  Taylor  series  expandions.  For  example,  in 
the  case  of  D{X),  write: 


The  effect  of  the  radius  and  mass  moment  of  inertia  of  the  rigid  disk  on  the  funda¬ 
mental  frequency  of  the  rigid  disk/flexible  beam  system  are  shown  in  figures  2  and  3. 
The  computation  of  the  exact  open-loop  transfer  functions  and  the  reduced  order 
open-loop  transfer  functions  using  a  Taylor  series  expansion  is  shown  in  appendix  A 
which  gives  a  Bode  plot  comparison  between  them.  The  beam  tip  response  and  disk 
angle  response  to  a  unit  step  torque  applied  to  the  disk  is  shown  in  figure  4.  Reduced- 
order  models  of  the  open-loop  transfer  functions  were  used  for  this  simulation  with  the 
first  four  vibration  modes  kept. 
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The  roots  of  N0(X)  are  the  same  as  the  roots  of  D(X)  for  e  =  0,  i.e.,  for  an  infinitely 
large  disk  moment  of  inertia.  Therefore,  the  colocated  zeros  [the  roots  of  N0(^)]  cor¬ 
respond  to  the  natural  vibration  frequencies  of  a  cantilevered  beam.  These  frequencies 
are  also  called  antresonance  frequencies.  If  the  actuator  is  driven  with  a  sinusoidal 
command  at  one  of  these  frequencies,  the  rigid  disk  will  not  move  and  the  elastic  beam 
vibrates  according  to  the  selected  cantilevered  mode  shape. 

The  tip  position  transfer  function  has  zeros  with  a  positive  real  part;  therefore,  the 
system  is  a  nnnminimum  phase  system.  The  initial  time  response  of  the  beam  tip  to  a 
unit  step  torque  input  to  the  disk  is  a  short  quasi-stationary  period  (approximately 
0.0066  seconds)  followed  by  motion  of  the  tip  in  the  opposite  direction  to  that  of  the  rigid 
disk  (fig.  4).  The  contributions  to  the  tip  position  of  the  rigid  body  mode  and  the  first 
flexible  modes.  This  type  of  initial  response  presents  a  difficult  control  problem. 

Reduced-Order  Model  with  Cantilevered  Modes 

It  is  common  to  use  the  constrained  mode  shapes  (cantilevered  modes)  rather 
than  the  unconstrained  mode  shapes  of  the  system  in  order  to  attempt  to  describe  the 
dynamic  behavior  of  an  elastic  structure.  Although  the  system  modes  lead  to  a  very 

simple  formulation  for  the  dynamic  equations,  they  are  difficult  to  use  for  a  multilink 

flexible  manipulator  because  of  their  dependence  on  the  relative  orientation  of  the  links. 
Formulation  in  tenns  of  the  cantilevered  modes  does  not  have  this  drawback. 

The  elastic  deflection  w(x,t)  is  expressed  in  terms  of  the  following  series: 

w(x,t)=^  <|)(x)  q(t) 

j  i  i 
i 

The  mode  shapes  for  a  cantilever  beam  have  been  previously  derived  in  various 
reference  books  and  textbooks  and  can  be  expressed  as  follows: 

0(x)  =  A  ■  [cos(fi  •  x)  +  cosh([3  ■  x)}  +  B  ■  [cos(|3  •  x)  -  cosh(|3  ■  x)  • 

+  C  •  [sin((3  •  x)  +  sinh(|3  x)  +  D  [sin([3  x)  -  sinh(|3  x)] 

The  boundary  conditions  are: 
at  x  =  0:  0(x)  =  0 

2  3 

at  x  =  L:  <D(x)  =  0  ^  q(x)  =  0  0(x)  =  0 

dx  dx 
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Therefore, 


A  =  C  -  O  and  g  = 


~COS([3 

L) 

+  cosh((3 

L)“ 

sin((3 

■  L)  -  sinh([3  •  L) 

sin(f3  ■ 

L) 

+  sinh([3  • 

L)_ 

cos  ((3 

L)  +  cosh(|3  •  L) 

Therefore,  the  mode  shapes  can  be  expressed  as: 


p(x)  =  -L 


cos((3  ■  x)  -  cosh((3  •  x)  - 


cosh  (A)  +  cos(A) 
sinh(A)TsTh(A)~ 


[sin(|3  ■  x)  -  sinh([3  ■  x] 


where  A  =  [3L  is  a  root  of  the  equation  [1  +  cos(A)cosh(A)]  =  0.  The  first  four  roots  of  this 
equation  are  1.8751,  4.6941,  7.8548,  and  10.9955,  and  the  values  of  D/B  correspond¬ 
ing  to  these  values  of  A  are  -0.7341,  -1.0185,  -0.9992,  and  -1.0000,  respectively.  The 
cantilevered  frequencies  are  related  to  the  values  of  A  by  the  following  equation: 


•'2  E  l  A  A 
or  =  — :  4  •  X 

0  L4 


In  the  expressions  for  kinetic  energy,  t,  and  potential  energy,  V,  in  w(x,t)  is  re¬ 
placed  with  its  series  representation.  The  result  is: 

T=  2  ■  [(i  +  Fi)  e2  +  e-£qi  |F2i  +  F3,|  +Sl4l2  F4,] 


v  =  5  Slqil2  F5i  -  u(t)  e 

1 


FI 


F2 


F3 


F4 


=  o  •  L  •  j^R2  +RL+- 

i  =  2  R  Jo  °  d(x)i  dx 
=  2  J*0  a  x  p(x)i  dx 

=  jo  ° '  !<>(x)*l2  dx 

4 


F5,=  LEI- 


2  0(x)i  dx 


]' 
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Lagrange’s  equations  was  applied: 


A  d 

dt  dq* 


d  T  d 
dqi  '  +  dqi 


V  =  0 


An  infinite  set  of  coupled  ordinary  differential  equations  was  obtained.  Using  the 
first  three  mode  shapes  of  the  cantilever  beam,  the  equations  of  motion  that  result  are 
as  follows: 


(j  +  F1)  B+l  £lF2i  +  F3i| q,  =u(t) 


\  •  |  F2,  +  F3, 1  •  0  +  F4|  ■  q’i  +  F5i  •  qf  =  0 


i  =  1  ,  2.  3 


There  are  four  equations  and  four  unknowns:  0  q1  q2  q3.  These  equations  can 
now  be  put  in  state-space  form  where  the  state  variables  are: 


x1  =  0 


x2  =  0 


X3  =  C1,  X5  =  Cl2 

x4  =  cii  x6  =  q2 

The  state-variable  equations  that  result  are  as  follows: 

i  =  1 , 2,  3,  4 


x7  =  q3 

x8  =  q3 


X2i-1=X2i 


X2  = 


1 _ 

J  +  FT 


1 


2  •  i  f  1 


J  +  FI 


F2i  +  F3j 
2  F4i 


[F2i  +  F3|| 
4  •  F4j 


U(t) 


-1 

'JTFT 


£  F5i  F2|  +  F3i 


r  F4i 


X2  i  +  1 


,  1  YTIF2,  +F3i 

1  "  J  +  FI  ^  4  •  F4j 

I  L_ 


.  .  F5j  F2,  +  F3i 
u(t)  +  I  — 5 - -  •  x2  ji  + , 


~F5~ 

F4| 

x2  i(1 


i  =  1. 2, 3 
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The  computation  of  the  stare-variable  equations  and  the  equations  of  motion  in 
state-space  matrix  form  are  shown  in  appendix  B.  Bode  plots  of  magnitude  (dB)  versus 
frequency  (rad/sec)  for  disk  angle  motion  and  beam  tip  motion  are  shown  in  figure  5. 
Also  shown  in  this  figure  are  the  beam  tip  response  and  disk  angle  response  to  a  unit 
step  torque  applied  to  the  disk. 

Finite  Element  Analysis 

A  finite  element  analysis  program  was  used  to  generate  a  dynamic  model  in  the 
following  form: 

[M]  [q]  +  [K]  [q]  =  [Q] 

where  [M]  is  the  mass  matrix,  [K]  is  the  stiffness  matrix,  [q]  is  the  nodal  displacement 
matrix,  and  [Q]  is  the  external  force  matrix.  This  model  is  valid  for  small  linear  elastic 
deformation  and  small  angular  velocities  of  the  rigid  base.  Using  the  same  parameter 
values  as  given  in  appendices  A  and  B,  the  output  from  the  finite  element  analysis 
program  is  given  in  appendix  C.  The  above  equation  can  be  transformed  into  a  set  of 
decoupled  modal  equations  as  follows.  First,  solve  the  eigenvalue  problem 


[M]  [U]  [co2]  =  [K]  [u] 

where 


[u]  =  modal  matrix 

[co2]  =  diagonal  matrix  of  eigenvalues 

The  response  may  be  described  as  a  superposition  of  normal  modes  in  the  form: 

[q]  =  [u]  m 

where  [p(t)]  is  a  column  matrix  consisting  of  a  set  of  time-dependent  generalized  coor¬ 
dinates.  Therefore, 

[q]  =  [u]  [Pi 

[M]  [up]  +  [K]  [u]  [p]  =  [Q] 

[u]T  [Ml  [u]  [PI  +  [u]T  |K]  [u]  [PI  =  [u]T  [O] 
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But  normal  modes  are  such  that 


[u]T  [M]  [u]  =  [I] 

[uf  [KJ  [u]  =  [0)21 

In  addition,  a  column  matrix  of  generalized  forces  N(t)0  was  introduced  associated 
with  the  generalized  coordinates  p(t)  and  related  to  the  forces  Q(t)  by 

[N]  =  [u]T  [Q] 


therefore, 


[fi  +  [«  ]  [PI  =  [N1 

This  represents  a  set  of  n  uncoupled  differential  equations  of  the  type 
P(t),  +  «  P(t),  -  N(t)r  r  =  1 . 2 . n 

Therefore,  modal  analysis  consists  of  uncoupling  the  equations  of  motion  by 
means  of  a  linear  coordinate  transformation;  the  transformation  matrix  is  just  the  modal 
matrix  [u].  Details  of  this  procedure  are  given  in  appendix  C.  Using  the  results  of  this 
analysis,  bode  plots  of  magnitude  (dB)  versus  frequency  (rad/sec)  for  disk  angle  motion 
and  beam  tip  response  and  the  beam  tip  response  and  disk  angle  response  to  a  unit 
step  torque  are  shown  in  figure  6. 

Lumped  Parameter  Model 

P.  Sheth  and  K.  Craig  have  developed  a  rigid  disk,  2  rigid  link,  2  torsional  spring 
model  of  the  rigid  disk/flexible  beam  dynamic  system  that  preserves  the  following 
characteristics  of  the  continuous  system: 

1 .  Natural  frequencies  of  the  first  two  modes  of  vibration 

2.  Mode  shapes  of  the  first  two  modes  of  vibration 

This  lumped  model  is  a  first  approximation  of  the  real  structure,  although  it  may  be 
an  excellent  approximation  for  specific  situations.  A  diagram  of  the  lumped  model  along 
with  its  parameters  is  shown  below. 
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Link  1:  ml,  LI,  II 


Rigid  disk:  R,  J 


Link  2:  m2,  L2,  12 


Parameter  values  (units:  in-lbf-sec): 

ml  =  0.000674  II  =  0.008346  kl  =  506.668  LI  =  18.96 

m2  =  0.000478  12  =  0.003814  k2  =  264.105  L2=  13.44 

a  =  3.6 

The  details  of  this  work  and  its  extension  to  other  multilink  dynamic  systems  are 
now  in  preparation  for  publication. 

Reduced-Order  Model  with  System  Modes 

A  standard  way  to  solve  the  dynamic  equations  previously  given  is  to  expand  the 
solution  y(x,t)  as  the  infinite  series: 


y(x,t)  =  Z  <!>(x)i  ■  q(t)i 
i 


where  <t>(x)  are  the  eigenfunctions  of  the  equation  of  motion  with  the  boundary  condi¬ 
tions,  and  q(t)  are  the  corresponding  time-dependent  generalized  coordinates.  0(x)  are 
called  the  system  or  unconstrained  mode  shapes.  To  find  the  mode  shapes  <Kx),  the 
external  torque  u(t)  is  set  equal  to  zero.  Then  search  for  the  eigensolutions  of  the 
equation  of  motion  and  boundary  conditions  of  the  form: 

j  w.  t 

i 

y(x,t).  =  <J)(x).  •  e 
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We  obtained  a  fourth  order  ordinary  differential  equation  in  the  variable  x  whose 
general  solution  is: 

<|>(x)  =  A  •  sin(p  x)  +  B  sinh((3  ■  x)  +  C  ■  cos(p  x)  +  D  cosh(P  •  x) 

This  equation  has  already  been  studied.  The  solution  to  the  following  system  of 
equation  is  given  in  appendix  D. 


-sin(X)  -  sinh(X)  -cos(X)  -  cosh(X)  R  •  cosh(X)  +  L  .  smh£X ) 

A 

-cos(X)  -  cosh(X)  sin(X)  -  sinh(X)  L  •  +  R  •  sinh(X) 

A 


2  ■  6  •  X 


-2 


L  e  X  -5  L  +  R 


A 

C 

0 


0 

0 

0 


D  =  R  0  -  C 


P  =  . 


-F 


£  = 


3 

G  L 


In  the  previously  given  expressions  for  kinetic  energy,  T,  and  potential  energy,  V, 
y(x,t)  =  (R  +  x)  0(t)  +  w(x,t)  with  its  series  representation  was  replaced.  The  result  is: 


-  i  ■  Z  4? 


J 


dx 


m 


|J  +  Jo  O-  <j>(x)fdx 


v  =  5  X  df  J  E  1  ^2  0(x)i  dx  -  U(t) 


dx2 


J 


■[?[ 


£  m 


■qi 


The 

d 

dx 


0(0) 


modes 

j-fo- 


were 

\2 


<t»(x)f  dx 


normalized  so  that  the 
equals  one  for  each  mode  shape. 


expression 


When  this  is  done,  the  expressions  for  T  and  V  reduce  to  the  following: 


T=i  Zq; 


V 


=  2  '  X  d?  ■  u(*)  ■  X  di 


adx  *i°)' 
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The  summation  is  over  the  range  0  to  3  where  the  value  i  =  0  corresponds  to  the 
rigid  body  mode  for  which  the  frequency  of  vibration  is  0.  Lagrange’s  equations  were 
then  applied,  and  a  set  of  decoupled  ordinary  differential  equations  of  the  following  form 
were  obtained: 

q,  +  q,  •  wf  =  u(t)  ~  ■  0(0),  i  =  o,  1 , 2,  3 

These  equations  are  easily  put  in  state-space  form  and  are  shown  in  appendix  E. 
The  Bode  plots  of  magnitude  (dB)  versus  frequency  (rad/sec)  for  disk  angle  motion  and 
beam  tip  response  and  the  beam  tip  response  and  disk  angle  response  to  a  unit  step 
torque  are  shown  in  figure  7.  These  results  are  slightly  different  from  previous  models 
as  slightly  different  values  of  the  parameters  J,R,L  and  d  were  used  to  match  actual 
experimental  parameters. 

Complete  Model  of  Electromechanical  System  and  Experimental  Verification 

The  reduced-order  model  with  system  modes  derived  in  the  previous  section  was 
used  as  the  model  for  the  mechanical  part  of  the  system.  The  following  equations  of 
motion  for  the  electrical  part  of  the  system  must  be  added. 

State  variables:  i,  and  i2 

i,  =  lj  [-R1  i,  -Kbl  0  +  Kal  -  e(t),] 

i2  =  -^2  [-R2  i2  -  Kb2  0  +  Ka2  •  e(t)2] 

The  torque  on  the  disk  from  the  motors  is  given  by: 

Torque  from  motors  =  Ktl  •  i,  +  Kt2  i2 

R1  =  electrical  resistance  of  motor  1 
R2  =  electrical  resistance  of  motor  2 
LI  =  electrical  inductance  of  motor  1 
L2  =  electrical  inductance  of  motor  2 

Kal  =  amplifier  gain  for  motor  1  (amplifier  1  was  linear  in  the  range  -1  <  0  <  1  input 
volts) 

Ka2  =  amplifier  gain  for  motor  2  (amplifier  2  was  linear  in  the  range  -1  <  0  <  1  input 
volts) 

Kbl  =  back  emf  constant  for  motor  1 
Kb2  =  back  emf  constant  for  motor  2 
e(t),  =  input  voltage  to  motor  1 

e(t)2  =  input  voltage  to  motor  2 
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To  this  torque  must  be  added  the  torque  due  to  Coulomb  friction,  Tp,  which  is 
represented  by  the  following  diagram: 


The  complete  equations  of  motion  for  the  electromechanical  system  are  given  in 
state-space  form  in  appendix  F  along  with  a  summary  of  all  the  experimentally  deter¬ 
mined  physical  parameters  of  the  dynamic  system. 

The  results  of  two  open-loop  simulations  for  the  complete  electromechanical 
system  along  with  an  example  of  the  block  diagram  used  in  MATRIXx  for  the  simulation 
are  shown  in  figure  8  through  1 1 . 

In  the  first  case,  the  input  from  both  motors  is  a  sine  wave  of  magnitude  0.5  volts 
and  frequency  8.5  Hz,  which  is  the  fundamental  frequency  of  the  mechanical  system. 
Only  the  output,  the  disk  angle,  from  this  simulation  is  shown. 

In  the  second  case,  the  input  from  both  motors  is  a  triangular  wave  of  magnitude 
0.125  volts  and  frequency  0.5  Hz.  Both  the  input  from  a  motor  and  the  output,  the  disk 
angle,  are  shown. 


CONCLUSIONS 

In  recent  years,  the  use  of  lightweight  materials  with  distributed  flexibility  in  ad¬ 
vanced  space  applications  and  in  the  construction  of  robotic  manipulators  has  led  to  a 
new  and  challenging  problem  in  the  combined  area  of  modelling  and  control.  Increas¬ 
ing  interest  has  arisen  to  properly  account  for  the  inherent  flexibility  of  these  structures, 
the  coupling  between  the  rigid  and  flexible  motions,  and  the  interaction  between  the 
control  system  and  the  flexible  structure  itself. 

In  this  report,  the  problem  of  modelling  a  flexible  beam  fixed  to  a  rigid  rotating  disk 
has  been  analyzed  in  great  detail.  The  equations  of  motion  and  boundary  conditions 
developed  by  a  means  of  Hamilton's  Principle  were  considered.  The  effect  of  the  radius 
and  mass  moment  of  inertia  of  the  rigid  disk  on  the  fundamental  frequency  of  the  sys¬ 
tem  was  observed.  Unconstrained  and  constrained  mode  expansions  were  obtained  for 
the  system.  In  the  absence  of  controller  dynamics,  it  was  found  that  the  constrained 
mode  expansion  yields  very  good  results  only  if  the  flexible/rigid  inertia  ratio  is  small. 
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Finite  element  and  lumped  parameter  models  of  the  system  were  also  derived.  The 
motor  dynamics  were  integrated  into  the  structural  model  considering  a  specific  labora¬ 
tory  setup.  This  experiment  was  also  used  to  validate  the  model. 

Multiple-link  robot  arms  with  motors  at  the  joints  can  be  considered  to  satisfy  the 
condition  of  small  flexible/rigid  inertia  ratio  so  that  constrained  expansions  are  expected 
to  perform  well.  However,  a  question  which  requires  further  research  effort  is  whether 
constrained  mode  expansions  together  with  controller  dynamics  are  suitable  for  cur¬ 
rently  projected  space  structures  that  have  very  large  flexible/rigid  inertia  ratios. 
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Figure  2.  Dynamic  analysis  of  a  flexible  beam  attached  to  a  rigid  rotor  with  variable  radius 
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TIP  RESPONSE  TO  UNIT  STEP  TORQUE  AT  T=0.005  SEC 


0  .  005  .  01  .015  .  02  .  025  .  03  .035  .04  .  045  .  05 


TINE  (SEC) 

DISK  ANGLE  RESPONSE  TO  UNIT  STEP  TORQUE  AT  T=0.005  SEC 


Figure  4.  Results  from  exact  analysis 
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Figure  5.  Results  from  reduced-order  model  with  cantilevered  modes 
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Figure  7.  Results  from  reduced-order  model  with  system  modes 
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APPENDIX  A 

RIGID  DISK  AND  FLEXIBLE  BEAM:  EXACT  OPEN-LOOP  TRANSFER 
FUNCTIONS  AND  REDUCED-ORDER  MODELS  OF  THE  OPEN-LOOP 

TRANSFER  FUNCTIONS 
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Parameter  values  (units:  inch-lbf-sec) 


O 

o 

CO 

II 

UJ 

L  :  =  36.0 

d  :  =  0.25 

8(R)  :  =  f 

p  :  =  0.000725 

a  d4 

64 

1  :  =  7t 

o  :  =  p  A 

E(J):=C  L= 

m :  -  Jr  ,  ■ 

6j(X)  :  =  2  ■  n  ■  f(X) 

R  :  =  5.625 

J  :  = 

0.416 

N(X)  :  =  1  +  cos(X)  •  cosh(X) 

NTIP(X,  R)  :  =  L  [sin(A.)  +  sinh(X)]  +  R  •  X  •  (cos(X)  +  cosh(X)] 


D(X,R,J)  :  =  e(J)  •  X3  •  [1  +  cos(X)  •  cosh(X)]  +  ( 1  +  X2  •  5(R)2|  sin(X)  ■  cosh(X) 


+  2  X  8(R)  •  sin(X)  ■  sinh(X)  +  |X2  ■  5(R)2  -  1 1  cos(X)  •  sinh(X) 


TF(X.R,J) 


-L  N(X) 

=  E  •  V  X  D(X,R,Jj 


TFTIP(X,R,J) 


-I  NTIP(X.R) 
E  •  I  X2  T5(X,FT,J) 


X  :  =  1.8  root[N(X),X]  =  1.875 


X  :  =  4.6  root[N(X),X]  =  4.694 


X  :  =  7.8  root[N(X),X)  =  7.855 


X  :  =  10.9  root[N(X),X]  =  10.996 


f(1 .875)  =  5.489 
(5(1 .875)  =  34.488 
f(4.694)  =  34.401 
co(4.694)  =  216.149 
f(7.855)  =  96.334 
63(7.855)  =  605.284 
f(  1 0.996)  =  188.78 
oi(10.996)  =  1.186  103 


X  :  =  2.3  root[D(X,R,J).X]  =  2.3993 


f(2.393)  =  8.941 
63(2.393)  =  56.176 
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X  :  =  4.8  root[D(X,R,J)A]  =  4.807 


f(4.807)  =  36.077 
63(4,807)  =  226.681 
X  :  =  7.8  root[D(X,R,J),X]  =  7.895  f(7.895)  =  97.31 8 

63(7.895)  =  611.464 

X  :  =  10.9  root[D(X,R, J) ,X]  =  11.018  f(1 1.018)  =  189.536 

63(1 1 .018)  =  1.191  ■  103 

X:  =  2.0 +  2.0  i 

root[NTIP(X,R),X]  =  2.12  +  2.1 2i  |2.12  +  2.12  i|  =  2.998 

X  :  =  5.0  +  5.0  i 

root[NTIP(\,R),>L]  =  5.081  +  5.081  i  |5.081  +  5.081  i|  =  7.186 

X.  :  =  8.0  +  8.0  i 

root[NTIP(X,R),X]  =  8.129  +  8.1 29i  |8. 1 29  +  8.129  i|  =  11.496 

X  :  =  11.0  +  11.0  i 

root[NTIP(X,R)A]  =  11.214  +  1 1 .21 4i  |1 1 .214  +  11.214  i|  =  15.859 


z  ^ 

^  II 

II  o 

PO 

63(2.998)  =  88.172 
63(7.186)  =  506.572 

N1(X):  =  ^N(>.) 

N1  {X)  =  0 

63(11.496)  =  1.296  103 

N2(X)  :  =  adxN1(X) 

N2(X)  =  -1 .91  ■  10  12 

63(15.859)  =  2.467  103 

N3(X)  :  =  ddx  N2(X) 

N3(X)  =  0 

D(A,R,J)  =  0 

D1  (X,R,J)  :  =  dd  D(X,R,J) 

D1(X,R,J)  =  1.563  •  10  11 

D2(X,R,J)  :  =  dd  D1(X,R,J) 

D2(X,R,J)  =  0 
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D3(X,R,J)  :  =  d°  D2(X,R,J) 
NTIP(X,R)  =  0 

NTIP1(A.,R)  :  =  ddxNTIP(X,R) 
N TIP2(X.R)  :  =  ^  NTIP1(X,R) 
NTIP3(X,R)  :  =  ^  NTIP2(X,R) 


NN(X)  :  = 


4  L2 


E  I  X2 


1  - 


DD(X)  :  -  o  327Q1T 


1  8744 

A4  ' 

"  2?3934 


D3(?i,R,J)  =  9.174 


NTIPI(X.R)  =  83.25 


NTIP2(A,R)  =  0 


NTIP3(>v,R)  =  -4.734  10 


1 1 


J.4 

4.6944 


1  - 


X4_ 

7.8554 


1  -  ^ 

10.9964 


1  - 


4.8074 


'l  .  _  *4J 

7.8954 


1  - 


X4 


11.0184 


NNTIP(^)  :  = 


166.5  L2 
E  I  •  X2 


1  + 


X4 


2.9984 


1  + 


X4 


7.1864 


1  + 


X4 

1 1 .4964 


1  + 


*4 


15.8594 
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APPENDIX  B 

REDUCED-ORDER  MODEL  OF  RIGID  DISK/FLEXIBLE 
BEAM  SYSTEM  WITH  CANTILEVERED  MODES 
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Parameter  values  (units:  inch-lbf-sec) 


E:  =  3.0  •  10 

2/  \  R 
5(r):  =  r 

e(J):  = 


a  O 


L:  =  36.0 


p:  =  0.000725 


d:  =  0.25 


I:  =  7t 
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A:  =  k 


F(X):  =  2: 


L 

■  k  V  O  L4 


X2 


ai(X):  =  2  k  1(X) 
R.  =  5.625 


J:  =  0.416 


Mode  shape  of  a  cantilever  beam:  <j>(x) 


XI 

X2 

X3 


=  1.8751 
=  4.69409 
=  7.85476 


Fa,  _  cosh(\)  +  cos(7i) 
'  F  ~  sjnh(X)  +  sin(Xy 


o(x,X)-  =  -L 


cos 


L  x 


-  cosh 


X 


Origin  =  1 


F(A.1)  =  0.734096 

F(X2)  =  1.018467 
F(X3)  =  0.999224 


■F(X) 


sin 


X 


-sinh 


x 


4>(L,\1)  =  71,999785 


§{L,X2)  =  -71.999917 


0(L,?i3)  =  72.000184 


F1:  =  oL 


R2  +  R  L  + 


L2' 


F3(X):  =  2  •  jo  x  •  <j)(x,  X)  ■  o  dx 


F2{X)\  =  2  •  R  ■  Jo  (J>(x,X)  •  a  dx 
F4(^):  =  Jo  o-[<D(x,X)]2  dx 


F5(X) 


4 


E  ■ 


d 

dx 


Sc 


-i2 


dx 


FI  =  0.853447 
F21  :=  F2(X.1) 


F21  =  0.406276 
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F31  :=  F3(X1 ) 
F41  :=  F4(X1) 


F31  =  1.888962 
F41  =  1.660399 


F51 

=  F5(M) 

F51  = 

1.975365  •  103 

F22 

=  F2(\2) 

F22  = 

0.225161 

F32 

=  F3(A.2) 

F32  = 

0.30142 

F42 

=  F4(X2) 

F42  = 

1 .660409 

F52 

=  F5(X2) 

F52  = 

7.758082  •  104 

F23 

=  F2(X3) 

F23  = 

0.132019 

F33 

=  F3(X3) 

F33  = 

0.107644 

F43 

=  F4(X3) 

F43  = 

1.660408 

F53 

=  F5(X3) 

F53  = 

6.082459  •  105 

State  variables 

:  xl  =  0 

x2  =  0dot  x3  =  ql 

x4  =  qldot 

x5  =  q2 

x6  =  q2dot  x7  =  q3 

x8  =  q3dot 

0  1 

0 

0  0  0 

0 

0  0 

3.206102 

103 

0  2.88881  •  104  0 

1.030814  105 

0  0 

0 

1  0  0 

0 

A  = 

0  0 

-3.405655 

■  103 

0  -1.99666  104  0 

-7.124681  104 

0  0 

0 

0  0  1 

0 

0  0 

-254.195034 

0  -4.901431  ■  104  0 

-8172785  103 

0  0 

0 

0  0  0 

0 

0  0 
k— 

-115.692072 

0  -1.042426  •  103  0 

-3.700429  •  105 

r 

ol 

!  2.348253  '' 

1 

i 

0 

r  i  ooo 

0  0 

|  -1.623043 

C  := 

41.625  0  72  0 

-72  0 

B  = 

0 

L 

,  -0.186181 

0 

ro^ 

L  -0.084737  j 

D  := 

lOJ 

APPENDIX  C 

FINITE  ELEMENT  ANALYSIS  OF  RIGID  DISK/FLEXIBLE  BEAM  SYSTEM 
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Using  six  beam  elements  to  represent  the  flexible  beam  offset  from  the  center  of 
rotation  of  the  disk  and  a  lumped  inertia  to  represent  the  rigid  disk,  the  mass  matrix,  [M], 
and  stiffness  matrix,  [K],  are: 


!  4.201  •  104 

14.94 

-8.729 

3.662 

-1.321  0.4059 

-0.1055 

1 

| 

17.52 

2.644 

-1.361 

0.5450  -0.1762 

0.04927 

18.52 

2.207 

-1.185  0.4371 

-0.1429 

18.70 

2.099  -1.009 

0.4293 

■  105  =  [M] 

18.88  1.662 

-1.227 

19.88 

3.345 

L 

5.992 

1.722  •  102 

-24.85 

9.041 

-2.413 

0.6413  -0.1599 

0.02659 

i 

j 

4.997 

-3.159 

1.271 

-0.3376  0.08415 

-0.01400 

\ 

3.895 

-2.863 

1.186  -0.2957 

0.04918 

i 

i 

3.811 

-2.821  1.102 

-0.1833 

o 

ro 

II 

25 

! 

3.727  -2.525 

0.6863 

! 

i 

2.625 

-0.9692 

; 

L 

0.4275 

J 

=  [Q] 

nodal  displacement  vector 


Torque 

0 

0 

0  =  [Q] 

0  external  force  vector 
0 

0  _ 


Solution  to  the  eigenvalue  problem:  [K]  [u]  =  [co2]  [M]  [u] 

C0q  =  0  6)5  =  3.1938  103  6J2  =  5.1292  104  w2  -  3.7616  105 

6)J  =  1 .4584  ■  1 06  w2  =  4. 1 546  •  1 06  co2  =  9.0099  •  1 06 
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2.469 

-2.535 

0.8549 

-0.3914 

0.2364 

-0.1620 

0.1389 

28.68 

-20.82 

-15.69 

46.52 

-67.04 

73.36 

-77.37 

43.34 

-12.27 

-54.78 

72.19 

-19.63 

-48.45 

89.98 

57.74 

7.786 

-71.04 

3.162 

69.07 

-0.9024 

-94.63 

71.91 

35.43 

-44.11 

-63.90 

-16.45 

50.69 

87.39 

85.97 

67.07 

20.09 

-21.95 

-52.03 

-64.88 

-66.79 

100 

100 

100 

100 

100 

100 

100 

Normal  modes  must  be  normalized  so  that:  [u]T[M]  [u]  =  [I] 

[j]t[K]  [u]  =  [oi2] 

The  result  is: 


0.8985 

-1.1314 

0.4669 

-0.2201 

0.1388 

-0.1002 

0.0660  " 

10.4357 

-9.2940 

-8.5689 

26.1614 

-39.3613 

45.3776 

-36.7462 

15.7682 

-5.4752 

-29.9182 

40.5986 

-11.5230 

-29.9682 

42.7346 

21.0086 

3.4755 

-38.7995 

1.7785 

40.5556 

-0.5582 

-44.9412 

26.1672 

15.8138 

-24.0894 

-35.9363 

-9.6564 

31.3545 

41.5060 

31.2821 

29.9402 

10.9742 

-12.3428 

-30.5504 

-40.1319 

-31.7223 

36.3838 

44.6376 

54.6172 

56.2393 

58.7158 

61.8603 

47.4935 

The  equations  of  motion  become:  [jj]  +  [652]  [(5]  =  [u]T[Q] 


0.8985 
-1.1314 
0.4669 
-0.2201 
0.1388 
-0.1002 
=  L0.066QJ 


[Torque] 


The  output  equations  are:  [q]  =  [u]  [[3] 

These  equations  are  uncoupled  and  it  is  straightforward  to  put  them  in  state  variable 
form. 
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APPENDIX  D 

REDUCED-ORDER  MODEL  OF  RIGID  DISK  -  FLEXIBLE  BEAM 
WITH  SYSTEM  MODES 
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NOTE:  Parameter  values  for  L,  d,  R,  and  J  are  different  here  than  in  previous  simula¬ 
tions.  These  are  the  actual  parameters  for  the  apparatus  on  which  experiments  were 
performed. 


E  :=  3.0  10 

L  :=  35.125 

S(R)  :=  ” 

p  :=  0.000725 

a  :=  p  AREA 

c(j)  o  L- 

co(X)  :=  2  •  n  •  f(X) 

J  :=  0.57301 

d  :=  0.249 

AREA  k  %  I  it  p 

FW  2~ n  » 


D1  (X,R,J)  :=  e(J)  •  X  3  •  [1  +  cos(X)  ■  cosh(X)]  +  1 1  +  X2  ■  5(R)2|  ■  sin(X)  ■  cosh(X).  .  . 


+  2  •  X  ■  8(R)  sin(X)  ■  sinh(X)  +  \X2  ■  5(X)2  -  1 )  cos(X)  sinh(X) 


D)X,R,J)  :  =  "2^  L  D1(X,R,J) 


X  :=  2.3 


X  :=  4.8 


X  :=  7.8 


X  :=  10.9 


XI  :=  root[d(X,R,J),X] 


X2  :=  root[D(X,R,J),X] 


X3  :=  root[D(X,R,J),X] 


X4  :=  root[d(X,R,J),X] 


XI  =  2.281132052 
f(X1)=  8.500007313 
63(X1 )  =  53.407121063 

X2  =  4.777449767 

f(X2)  =  37.282987411 
w(X2)  =  234.255918709 

X3  =  7.884939741 

f(X3)  =  101.558248674 
co(X3)  =  638.109295891 

X4  =  11.012274831 
f(X4)  =  198.094550571 
65(X4)  =  1.24466477  103 


X1 1  (X)  :=  -sin(X)  -  sinh(X) 

X12(X)  :=  -cos(X)  -  cosh(X) 

XI  3(X)  :=  R  ■  cosh(X)  +  -  '  Sl^h(X) 


X21  (X)  := 
X22(X)  := 

X23(X)  := 


-cos(X)  -  cosh(X) 
sin(X)  -  sinh(X) 

L  •  cosh(X) 


X 


+  R  •  sinh(X) 
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X31(X,R)  :=  2  8(R)  •  X 
X32  :=  -2 

X33(X,J,R)  :=  -  L  •  8(R)  +  R 


'  XI  1(A.)  X12(X) 

X(X,R,J)  :=  X21  (X)  X22(X) 

_X31(X,R)  X32 

X13(X)~ 

X23(X) 

X33(X,J,R)_ 

X  :=  XI 

X  = 

2.281132052 

7-5.600935951 
X(X,R,J)  =  -4.2928731 

u  0.787437043 

-4.2928731 

-4.084649959 

-2 

104.548427761" 

105.502227379 

68.454784401. 

|X(X,R,J)|  =  1.370693935  • 

10  5 

X  :=  X2 

X  = 

4.777449767 

7-58.398403556 
x(X,R,J)  =  -59.469720184 

L  1.649155258 

-59.469720184 

-60.394172143 

-2 

796.837320401" 
796.848176783  j 
300.258431 662 j 

|x(X,R,J)|  =  -6.553719071  ■ 

10  8 

X  :=  X3  X  =  7.884939741 


1^1 .329  481 791  1 03  -1 .328451 693  •  103  1.397190909  104 

X(X,R,J)  =j  -1 .328451 693  1 03  -1.327482749  103  1.397190849  104 

|_  2.721847526  -2  817.899062996  J 

|X(X,R,J)|  =  -4.009991272  107 

The  solution  to  the  set  of  homogeneous  equation  with  0.  =  1  is  given  by  the  following: 

ORIGIN  =  1 

i  :=  1,2.  .3  rxr  r  31.938591] 

X  :=  X2  .  C  :=  j  87.720856  i 

_X3j  79.353947 J 


46 


r  -5.81 3294^ 

A  :=  I  -75.685245 
L- 168. 705692  j 

The  modes  of  vibration  are: 
<j)1(x)  :=  A,  ■  sin 


02(x)  :=  A2  •  sin 

03(x)  :=  A3  ■  sinj 
0O(x)  :=  R  +  x 


"  21.211351 
B  :=  83.037494 

Li  73. 160387. 


r  -25.876091 
D  :=  :  -81.658356, 
L-1 73.291 447J 


~x, 

+  Bi  ■  sinh 

r*i  vi 
r  x 

r .  -i 

X, 

L”  '  x 

+  C,  ■  cos 

L  x 

+  Di  cosh 

L  x 

r~  .  ~~i 

h  .  x 

L  x 

+  B2  sinh 

i  r 

X 

_ 1  L 

+  C2  ■  cos 

h  .  x 

L  x 

+  D2  cosh 

i  r 

X 

^!_j 

_ 1  L 

^3  .  x 

L  x 

+  B3  ■  sinh 

1 

X 

001  | 

I-1 

_ 1 

+  C3  ■  cos 

h  . 

u  L 

+  D3  cosh 

1 

X 

_ i 

This  is  the  RIGID  BODY  mode. 
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APPENDIX  E 

STATE-SPACE  FORMULATION  OF  EQUATIONS  OF  MOTION  FOR 
REDUCED-ORDER  MODEL  OF  RIGID  DISK  -  FLEXIBL  E  BEAM 
SYSTEM  WITH  SYSTEM  MODES 


x  =  Ax 

+  Bu 

y  = 

Cx  + 

Du 

To 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

-2 

-CO 

0 

0 

0 

0 

0 

=  [A] 

0 

0 

0 

0 

o2 

1 

0 

0 

0 

0 

0 

0 

-C02 

0 

0 

0 

0 

0 

0 

0 

0 

0 

02 

1 

0 

0 

0 

0 

0 

0 

-"3 

0 

OH 

(ML) 

1 

<ML) 

2 

<ML)3 

0 

0 

0 

£0(°)0 

0 

0 

£0(0)2 

o  4  0(0)3 

0 

Si 

0 

ft  0(0),  =  IB] 

0 


a;  0(0)2 
0 


d 

dx 


<>(0)3  i 

-J 


1 


'O' 

.0. 


=  [D] 


Output  1 :  absolute  beam  tip  position  (inches) 
Output  2:  disk  angle  (radians) 

Input:  Torque  u(t)  to  disk 
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APPENDIX  F 

COMPLETE  EQUATIONS  OF  MOTION  IN  STATE-SPACE  FORM  FOR  THE 

ELECTROMECHANICAL  SYSTEM 


oo  o  o  oooo 


State  of  variable  equations:  x  =  Ax  +  Bu 


y  =  Cx  +  Du 


0 

0 

1 

-2  •  T1  0), 
0 


=  columns  1  through  5  of  [A] 


-Kbl  d  x/fn  A 

Tr  ®  '*  0 


-Kb2  d  ,n.  n 
~L2~  dx  <t)(0)o  0 


-Kbl  d  ,/m  n 

-XI”  c5f M°h  0 

-Kb2  d  n 

TT'di**01,  0 


Kll  ddx-  0(0),  Kt2  ^  <|)(0)o 


Ktl  ^flO),  Kt2  ^  <1.(0), 


-2n,o>, 


Ktl  ^  *(0)2  Kt2^  4(0)2  • 


=  columns  6 
through  10 

of  [A] 


-co3  -2Ti3a)3 


Kt1^  <M0)3  Kt2  ~<t>(0)3 


-Kbl  d 


-Kbl  d 


tt-cE* 0)2  0  tv 

^■^MOh  o  f  |#)3  o 
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To 

0 

01 

:  0 

0 

o!  =  [D] 

LO 

0 

OJ 

Inputs: 

1  -  el  (t)  =  voltage  input  to  motor  1 

2  -  e2(t)  =  voltage  input  to  motor  2 

3  -  Tp  =  Coulomb  friction  torque  on  disk 

Outputs:  1  -  beam  tip  position  (inches) 

2  -  disk  angular  position  (rad) 

3  -  disk  angular  velocity  (rad/sec) 

Summary  of  parameter  values:  (units:  in-lbf-sec) 


E  :=  3.0  10' 

L  :=  35.125 

p  :=  0.000725 

R  :=  6.0625 

J  :=  0.573 

d  :=  0.249 
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R1  :=  29  LI 


R  :=  55  Q 


LI  :=  0.0125  H 


L2  :=  0.025  H  Ka1:=  17.84  Ka2:=  17.92 

Ktl  :=  4.204  in-lbf/A  Kt2  :=  5.51 1  in-lbf/A 

Kbl  :=  0.473  V-sec/rad  Kb2  :=  0.620  V-sec/rad 

n,  :=  .01  n2  :=  01  n3  :=  01 

Tg  0.39  in-lbf 


E  =  modulus  ot  elasticity  of  beam  material 

L  =  length  of  beam 

p  =  mass  density  of  beam  material 

R  =  radius  of  rigid  disk 

J  mass  moment  of  inertia  of  rigid  disk 

d  =  diameter  of  uniform  beam 

R1  =  electrical  resistiance  of  motor  1 

R2  =  electrical  resistance  of  motor  2 

LI  =  electrical  inductance  of  motor  1 

L2  =  electrical  inductance  of  motor  2 

Kal  =  amplifier  gain  for  motor  1 

Ka2  =  amplifier  gain  for  motor  2 

Ktl  =  motor  1  torque  constant 

Kt2  =  motor  2  torque  constant 

Kbl  =  motor  1  back  emf  constant 

Kb2  =  motor  2  back  emf  constant 

q  =  damping  coefficient  associated  with  each  mode 

Tp  =  Coulomb  friction  torque  acting  on  disk 

oi  =frequency  in  rad/sec 


<t>(L)0  :=  34.9017 
ddx-<|)(0)o  :=  0.8474 
co,  :=  53.4071 


(ML),  :=  -46.5615 
^ -0(0)1  :=  0.9226 
c'd2  :=  234.2559 


0(L)2  :=  55.7653 
^■0(0)2  :=  0.3433 
di  :=  638.1093 


0(L)3  :=  -56.5945 
g* -0(0)3  :=  0.1625 
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